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PREFACE 


The  use  of  underground  cavities  to  muffle  or  seismlcally  de¬ 
couple  nuclear  explosions  has  been  the  subject  of  considerable 
theoretical  effort  at  RAND.  Previous  work  has  suggested  that  the 
use  of  such  decoupling  can  reduce  the  transmission  of  seismic  sig¬ 
nals  from  underground  nuclear  detonations  by  factors  of  several 
hundred,  depending  on  such  factors  as  the  nature  of  the  cavity  wall 
pressures  thus  generated.  This  memorandum  reports  on  the  wall  pres¬ 
sures  and  temperatures  derived  from  theoretical  calculations  of  the 
radiation  and  hydrodynamic  actions  from  1.7  KT  explosions  in  20  and 
40  meter  cavities  filled  with  air  at  various  densities. 

The  calculations  were  made  in  support  of  the  testimony  of 
A.L.  Latter  at  the  Geneva  Nuclear  Test  Ban  Conference,  December  1959, 
and  the  theory  developed  by  Latter,  LeLevier,  Martinelli,  and  McMillan 
as  reported  in  1961. 

Interest  in  decoupling  calculations  lagged  somewhat  during  the 
active  nuclear  testing  period  (1961-1963)  that  followed  the  three 
year  moratorium  (1958-1961).  The  signing  of  the  partial  test  ban 
treaty  and  the  further  attempts  to  extend  coverage  to  a  total  test 
ban  hive  provided  renewed  interest  in  technical  details  of  cavity 
explosions  such  as  provided  by  these  numerical  results.  In  addition, 
the  results  have  found  increasing  application  to  possibilities  for 
underground  testing  where  large  cavities  might  be  used  to  explore 
close-in  fireball  effects  and  high  temperature,  high  pressure  physics. 
Recently,  several  organizations  engaged  in  theoretical  planning  for 
future  underground  tests  have  expressed  interest  in  these  calcula¬ 
tions.  It  is  hoped  that  this  Memorandum  will  be  directly  applicable 
to  these  activities. 
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SUMMARY 


Underground  cavities  have  been  suggested  as  an  effective  means 
of  reducing  the  seismic  signals  from  nuclear  detonations.  This  memo¬ 
randum  reports  the  wall  pressures,  temperatures,  and  densities  for 
five  specific  examples,  for  which  average  cavity  pressures  range 
between  44  and  420  bars,  with  peak  pressures  between  0.2  and  8.2  Kbars. 
One  nuclear  yield,  1.7  kilotons,  was  exploded  (theoretically)  in  two 
cavities  of  20  and  40  meters  radius,  each  with  normal  and  one- tenth 
normal  air  density  in  it.  The  larger  cavity  was  also  used  with  a 
density  of  one-hundredth  normal.  In  addition  to  the  wall  pressures, 
temperatures,  and  densities  as  functions  of  time  for  each  case,  some 
typical  pressure,  -  temperature,  and  density  profiles  in  the  cavity 
interiors  at  various  times  are  reported. 


T 


-vii- 


CONTENTS 

PREFACE . iii 

SUMMARY .  v 

Section 

I.  INTRODUCTION .  1 

II.  PHYSICAL  AND  MATHEMATICAL  MODEL .  2 

III.  RESULTS .  5 

REFERENCES .  11 

FIGURES . 12 


-1- 


I.  INTRODUCTION 


Following  the  early  theoretical  work  on  seismic  decoupling  of 

(1-4) 

underground  nuclear  explosions,  experimental  efforts  largely 

confirmed  the  suggested  range  of  decoupling  factors.^  Such  infor¬ 
mation,  combinec  with  the  data  on  transmission  of  seismic  signals  from 
explosive  sources,^  ^  has  been  quite  germane  to  the  questions  of 
inspection  and  enforceability  of  a  more  complete  test  ban  treaty  and 
the  detection  of  clandestine  tests.  The  theoretical  calculations 
reported  or>  in  this  memorandum  were  undertaken  in  1958  as  a  part  of 
further  investigation  of  the  decoupling  phenomena  suggested  and 
reported  in  Ref.  1.  These  calculations  define  more  exactly  the  pres¬ 
sures  and  temperatures  to  which  a  cavity  wall  might  be  exposed  for 
one  yield  (1.7  KT) ,  two  cavity  radii  (20  and  40  m) ,  and  three  different 

if 

air  densities  (normal,  one- tenth  normal  and  one  hundredth  normal). 

The  physical  assumptions  and  limitations  which  have  governed 
these  numerical  calculations  are  briefly  outlined  and  justified  in 
Section  II,  along  with  numerical  procedures  and  the  mathematical 
model.  The  graphical  results  are  presented  in  Section  III  together 
with  some  interpretive  remarks. 


The  20  m/ one- hundredth  normal  density  calculation  was  not  com¬ 
pleted  and  is  not  reported  because  the  very  small  anount  of  air  in 
that  sphere  did  not  significantly  alter  the  wall  pressures  from  those 
expected  in  a  wholly  evacuated  sphere. 
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II.  PHYSICAL  AND  MATHEMATICAL  MODEL 

The  basic  physical  assumptions  of  spherical  symmetry,  thermo¬ 
dynamic  equilibrium,  radiation  diffusion,  and  hydrodynamics  have  been 
shown  to  be  quite  valid  for  the  early  phases  of  a  nuclear  explosion 
in  the  atmosphere.  ^  At  late  times,  when  gravity  has  had  time 
to  distort  the  flow  or  at  times  when  the  motions  have  slowed  consider¬ 
ably  allowing  instabilities  to  grow,  the  spherical  symmetry  may  not 
be  maintained.  The  problems  of  interest  here  are  not  subject  to  such 
distortions  until  well  after  the  shock  reaches  the  cavity  walls,  but 
the  late  time  reverberations  reported  in  Section  IV  are  likely  un¬ 
realistic  to  the  extent  that  instabilities  and  turbulence  have  des¬ 
troyed  some  of  the  symmetry.  Gravity,  however,  could  not  have  con¬ 
tributed  much  in  times  less  than  a  second. 

A  further  assumption  which  requires  some  examination  is  that  of 
no  wall  response.  The  walls  are  here  taken  to  be  rigid  reflecting 
surfaces  and  with  very  high  opacity,  so  that  no  shock  energy  and  no 
radiation  flows  into  them.  If  the  walls  are  hit  hard  enough  they 
must  move  and  become  heated  (and  even  at  low  shock  inputs  the  walls 
will  transmit  some  elastic  strain  energ$,  so  that  whether  the  walls 
respond  or  not  is  not  the  question,  but  rather  whether  the  response 
has  any  significant  influence  on  the  phenomena  inside  the  cavity. 

To  influence  the  hydrodynamics  appreciably,  the  walls  would  have  to 
move  in  such  a  way  that  the  shape  or  volume  of  the  cavity  was  changed 
by  a  few  per  cent.  The  consequences  of  lesser  changes  in  radius  or 
shape  would  be  difficult  to  distinguish  in  the  results  of  such  a 
finite  difference  numerical  solution.  In  a  rock  cavity  large  enough 
that  the  walls  remain  elastic  in  responding  to  contained  explosion, 
the  wall  motions  are  negligible. 

Even  while  the  gross  response  of  the  walls  may  be  elastic  and 
involve  only  small  strains,  there  remains  some  question  of  the  in¬ 
fluence  of  cracks,  spalls,  rock  falls,  and  wall  melt  or  vaporization. 
Spalling  seems  least  likely  to  occur  early  enough  to  affect  the  first 
reflected  pulse  or  two.  Except  for  rather  complicated  geological 
features  in  the  cavity  vicinity  and  for  odd  shaped  cavities  or  abnormal 
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shock  loads,  mechanisms  for  spall  generation  do  not  seem  to  be  present. 

Cracks  may  propagate  into  the  rock  and,  on  opening  up  under  the 
stress  in  the  cavity,  accept  an  appreciable  flow  and  volume  of  cavity 
gas.  In  such  an  event,  the  first  pulse  reflection  is  not  likely  to 
be  influenced,  but  subsequent  cavity  pressures  and  temperatures  could 
be  greatly  reduced. 

To  open  up  an  appreciable  volume  of  crack  space  would,  however, 
require  that  cavity  pressures  at  relatively  late  times  remain  appre¬ 
ciably  higher  than  either  the  strength  (in  tension  or  shear)  of  the 
rock  or  the  overburden  pressure.  Cracks  will  not  torro  if  the  walls 
remain  elastic  and,  even  if  the  walls  flow  plastically,  cracks  cannot 
open  until  the  hoop  stress  exceeds  the  compressive  stress  of  the 
expanding  explosive  waves.  Hoop  stresses  are  developed  only  as  the 
walls  expand  radially  and  can  propagate  cracks  only  in  regions  where 
radial  motions  have  been  sufficient  to  raise  those  hoop  stresses 
above  the  compressive  stress  plus  that  of  the  overburden  pressure. 

Cavity  collapse  due  to  rock  falls  and  chimney  formation  is  always 
a  slow  process  compared  to  the  shock  propagation  and  reflection  times, 
and  no  account  seems  necessary  of  the  effect  of  such  phenomena  on  the 
cavity  wall  pressures  and  temperatures  within  the  first  second. 

Wall  melt  and  vaporization  co-  Id  cause  the  most  serious  deviation 
of  the  cavity  phenomena  from  that  predicted  on  the  basis  of  no  wall 
interaction.  Vaporizing  even  a  rather  thin  layer  of  wall  material 
can  alter  appreciably  the  amount,  temperature,  or  pressure  of  the 
cavity  gas.  It  can  also  appreciably  alter  the  opacities  and  so  the 
radiative  energy  transfer  properties  of  the  cavity  vapors.  As  sug¬ 
gested  in  Ref.  1,  pressure  due  to  vaporization  at  the  wall  is  entirely 
negligible  below  one  electron  volt,  and  may  not  be  serious  at  even 
12  ev.  In  these  examples,  wall  temperatures  fall  within  or  below 

this  range.  Fragmentary  evidence  such  as  that  from  the  Cowboy  Program 

(12) 

experiments  shows  good  agreement  between  calculations  and  data  '  7 

and  indicates  that  wall  interactions  do  not  seriously  affect  the 

cavity  phenomena  for  the  first  several  reverberations.  The  Cowboy 

(12) 

comparisons  were  generally  at  lower  peak  pressures  and  from  high 
explosive  sources,  however. 
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With  the  exception  of  the  rigid  boundary  which  represented  the 
cavity  wall  and  the  omission  of  the  thermal  radiation  loss  consider¬ 
ations,  these  calculations  used  the  same  numerical  procedures  and  the 
same  machine  program  employed  for  several  free  air  nuclear  burst  prob¬ 
lems.  The  model  is  best  described  in  Ref.  10,  in  which  the  expression 
of  the  hydrodynamic  and  radiation  diffusion  mechanisms  as  differential 
equations  in  Lagrangian  form  are  expressed  in  finite  difference  form, 
and  procedures  for  step-wise  integration  are  outlined. 

Equations  of  state  and  opacities  were  considered  in  some  detail 
not  only  for  the  air  in  the  cavities  but  for  the  various  bomb  vapors 
as  well.  The  equation  of  state  for  air  was  included  in  an  analytic 
form  which  approximated  the  theoretical  thermodynamic  data  to  within 
107.  at  almost  all  temperatures  and  densities. 

The  weapon  was  assumed  to  be  composed  of  spherical  shells  of 
appropriate  vapors  in  the  interior  of  which  the  1.7  KT  of  energy  was 
suddenly  released.  Since,  even  in  the  lowest  density  cases,  the  bomb 
vapors  do  not  expand  to  much  more  than  one  half  the  cavity  radius,  it 
is  the  air  rather  than  the  bomb  material  which  is  of  most  concern  in 
describing  cavity  wall  loads.  However,  mixing  and  jetting  can  eventually 
bring  the  hot  bomb  vapors  in  contact  with  the  walls,  where  the  conden¬ 
sable  matter  would  in  the  end  be  deposited. 
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III.  RESULTS 


Five  separate  problems  are  reported,  all  for  the  1 . 7  KT  weapon 
in  either  20  meter  or  40  meter  radius  cavities.  Air  densities  of 
normal  (p^  =  1.293  gm/ liter)  and  one- tenth  and  one-hundredth  normal 
were  used.  The  following  table  summarizes  some  of  the  parameters 
observed  at  the  walls  for  the  various  cases. 

Peak  pressures  on  the  wall  do  not  follow  a  simple  scheme. 
Although  doubling  the  radius  drops  the  wall  pressure  by  eleven- fold 
at  normal  air  density,  it  drops  by  only  a  little  more  than  a  factor 
of  three  for  the  one- tenth  normal  density  case.  For  a  fixed  cavity 
size,  the  decrease  of  ambient  density  and  pressure  does  not  guarantee 
a  drop  in  peak  pressure.  The  average  pressure  on  the  wall  shows 
hardly  any  dependence  on  the  ambient  initial  density. 

The  first  figure  indicates  the  regularity  with  which  the  rever¬ 
berations  progress  in  the  normal  density  20  meter  cavity  case.  Wall 
pressure  as  a  function  of  time  is  shown  in  Fig.  2.  At  the  time 
labeled  "doubled  zones,"  the  running  time  of  the  numerical  problem 
was  increased  by  coarsening  the  mass  zone  system  by  a  factor  of  two. 
Such  a  sudden  change  in  the  mass  sizes  introduces  some  disturbance, 
and  the  subsequent  history  shows  some  loss  in  detail.  Wall  tempera¬ 
tures  for  this  case  are  shown  in  Fig.  3,  and  air  density  at  the  wall 
is  included  in  Fig.  4.  Pressure,  temperature,  and  density  all  show 
the  same  time  history  details,  and  the  wall  conditions  are  clearly 
dominated  in  this  case  by  shocks  and  hydrodynamic  motions. 

Pressure  profiles  for  this  first  case  are  shown  for  times  before 
the  first  reflection  in  Fig.  5.  The  shock  formation  occurs  at  around 

6  meters,  and  is  more  clearly  illustrate^  r.  subsequent  figure 
(Fig.  12).  The  complications  that  arise  from  secondary  shocks  is 
evident  in  this  figure  -  Fig.  5.  A  case  shock  representing  the  energy 
of  the  bomb  debris  overtakes  the  main  shock  just  as  it  reaches  the 
wall  and  further  complicates  the  reflection  there.  The  reflection 
factor  of  around  twelve  is  to  be  expected  at  that  incident  pressure 
(700  bars)  . 

After  reflection  the  interior  zones  are  so  expanded  that  the 
shock  no  longer  has  a  sharp  front  in  a  radial  profile  plot  (Fig.  6)  , 


Table 


6 


X? 

CU 

D 

_c 

■+- 

, 

O 

u 

a 

<u 

CL 

E 

3 

QJ 

4- 

4- 

o 

l_ 

_ 

4> 

wal 

Q. 

E 

<u 

05 

C 

c 

tn 

0 

o 

L. 

4— 

u 

c 

<4- 

o 

<y 

4- 

u 

JZ 

0) 

•t 

u- 

3 

0) 

t2 

c 

u 

1 

o 

4- 

LO 

1 

e 

o 

UO 

•4— 

o 

"a 

1 


After  many  reflection,  the  wall  temperature  rises  toward  an  asymptotic  value  of  nearly  three  electron 
volts  (See  Fig.  35). 

Although  the  cavity  rapidly  approaches  a  nearly  isothermal  state,  the  temperature  seems  to  oscillate  more 
than  twice  as  rapidly  as  the  pressure  and  density  do  -  a  consequence  of  the  radiation  dominance  in  this  case. 
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but  the  imploding  shock  wave  can  be  seen  to  grow  in  strength  as  it 
again  approaches  the  origin.  The  shock  is  illustrated  at  times  after 
the  first  reflection  at  the  origin  but  before  the  second  wall  reflec¬ 
tion  in  Fig.  7.  At  subsequent  times,  the  interior  of  the  cavity  grows 
more  nearly  uniform  in  pressure,  with  relatively  weak  shocks  batting 
back  and  forth. 

Temperature  profiles  for  the  first  case  are  shown  prior  to  first 
reflection  in  Fig.  8  and  at  late  times  in  Fig.  9.  The  high  tempera¬ 
ture  gradient  beyond  15  meters  would  undoubtedly  not  hold  up  for  long 
under  the  influence  of  instabilities  and  mixing.  Early  and  late  den¬ 
sity  profiles  are  similarly  illustrated  in  Figs.  10  and  11. 

The  early  explosion  history  is  dominated  by  radiation  diffusion 
ahead  cf  the  bomb  shock.  For  some  distances  the  diffusion  growth  is 
faster  than  a  shock  wave  at  the  same  air  temperature,  but  the  diffusion 
growth  slows  as  the  driving  temperature  drops,  and  the  rate  of  feeding 
energy  into  the  radiation  wave  decreases.  A  shock  wave  gradually 
builds  up  at  this  front,  and  eventually  dominates  the  expansion. 

Figure  12  illustrates  this  transition,  showing  the  gradual  increase 
of  overdensity  at  the  front.  Times  and  radii  are  indicated  for  both 
the  overdensity  and  the  front  temperature.  Also  shown  is  the  interior 
temperature  at  corresponding  times  (but  not  as  a  function  of  the 
radius)  . 

For  the  case  of  a  40  m  cavity  at  normal  density,  the  shock  posi¬ 
tion  versus  time  is  shown  in  Fig.  13;  and  pressures,  temperatures, 
and  densities  on  the  wall  are  given  as  functions  of  time  in  Figs.  14, 
15,  and  16,  respectively.  Again,  the  detailed  shock  structure  is 
evident  in  these  time  histories  at  the  wall,  and  are  more  obviously 
related  to  multiple  shocks  in  the  pressure  profiles  shown  in  Figs.  17, 
18,  and  19.  Corresponding  temperature  profiles  (Figs.  20  and  21)  and 
density  profiles  (Figs.  22  and  23)  continue  to  demonstrate  the  basic 
hydrodynamic  nature  of  such  a  large  cavity  filled  with  norma]  air. 

Pressure,  temperature,  and  density  at  the  wall  are  shown  as 
functions  of  time  for  the  20  m  cavity  evacuated  to  one- tenth  normal 
in  Figs.  24,25,  and  26.  The  wall  pressures  quickly  danp  towards  an 
average  value  (Fig.  24)  of  slightly  more  than  400  atmospheres. 
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The  damping  of  the  reverberations  may  be  in  part  due  to  the  con¬ 
tinued  radiation  flow  demonstrated  in  Fig.  25.  The  rather  rapid 
rise  of  wall  temperature  after  the  first  reflection  is  a  consequence 
of  this  continued  heating.  Pressure  profiles  show  the  same  influence 
(Fig.  27)  in  the  slow  build-up  of  a  shock  and  the  rather  low  reflec¬ 
tion  factor  (about  3.3).  Late  time  pressure  profiles  show  relatively 
small  pressure  fluctuations  (Fig.  28) .  The  temperature  profiles  for 
this  case  (Figs.  29  and  30)  show  the  initial  radiative  dominance 
(Fig.  29),  and  the  modest  changes  caused  by  shocks  at  late  times 
(Fig.  30)  .  Density  profiles  (Fig.  31)  at  early  times  appear  to  be 
quite  complex,  because  of  the  slow  development  of  the  air  shock,  and 
the  more  dominant  excursion  of  the  bomb  or  case  shock.  The  late  time 
densities,  however,  become  rather  uniform  (Fig.  32). 

For  the  larger  (40  m)  cavity  at  one- tenth  normal  air  density, 

Fig.  33  shows  a  plot  of  shock  positions  as  a  function  of  time,  while 
Fig.  34  describes  the  wall  pressure  versus  time  for  many  reverbera¬ 
tions.  With  less  air  in  the  cavity,  the  radiation  continues  to  dif¬ 
fuse,  and  Fig.  35  shows  its  influence  superimposed  on  the  shock  ef¬ 
fects  on  wall  temperature.  A  corresponding  but  opposite  drift  must 
exist  in  the  air  densities  near  the  wall,  since  the  pressures  at  the 
wall  exhibited  no  obvious  change  from  the  inter-cycle  /erage,  and  it 
does  show  in  Fig.  36.  The  pressure  profiles  for  this  case  (40m, 
l/10th  normal)  show  some  of  the  irregularities  introduced  by  secondary 
shocks  as  well  as  the  progress  of  the  main  air  shock  and  its  reflec¬ 
tions  (Figs.  37  and  38). 

The  temperature  profiles  of  Figs.  39  and  40  show  features  similar 
to  the  corresponding  plots  for  the  normal  density  case  (Figs.  20  and 
21),  but  with  smaller  temperature  gradients  sustained  for  shorter 
times.  The  density  pi  files  (Figs.  41  and  42)  show  the  shock  forming 
at  a  larger  radius  than  in  the  normal  density  case  (Figs.  22  and  23) 
but,  in  keeping  with  the  temperature  differences,  developing  a  shallow¬ 
er  fireoall.  In  this  case,  the  interior  densities  are  mostly  greater 
thas  one-fifth  of  the  initial  ambient  cavity  density,  while  in  the 
previous  example  with  normal  density  air,  the  fireball  interior  drops 
to  less  than  one-tenth  of  the  initial  ambient  value. 
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When  the  40  m  cavity  is  pumped  down  to  one-hundredth  normal,  the 
phenomena  at  the  wall  changes  character.  The  wall  first  sees  the 
radiation  diffusion  wave,  and  the  wall  pressure  (Fig.  43)  as  well  as 
temperature  (Fig.  44)  show  a  slow  rise  over  the  first  half  millisecond 
before  the  shock  strikes  the  wall.  Actually,  the  radiation  front  is 
slowing  down  and  beginning  to  shock  up,  so  that  the  density  at  the 
wall  shows  some  slight  rise  before  shock  arrival  (Fig.  45).  Perhaps 
more  interesting,  however,  is  the  extra  oscillation  that  persists  in 
the  temperature  profiles.  Every  other  pulse  in  Fig.  44  (wall  tempera¬ 
ture  versus  time)  corresponds  to  the  shock  reflection,  but  the  com¬ 
parable  temperature  pulses  in  between  have  apparently  no  accompanying 
consequence  in  either  the  pressure  or  the  density.  Of  course,  there 
must  be  some  consequence,  but  the  temperature  variations  are  not 
large  to  begin  with,  and  the  shocks  are  nearly  isothermal. 

The  pressure  profiles  for  this  case  show  the  diffusion  front 
preceding  the  shock  wave  all  the  way  to  the  wall  (Fig.  46).  The 
subsequent  reflected  shocks  are  quite  weak  (Fig.  47)  except  as  they 
implode  on  the  origin  each  time.  The  temperature  profiles  shown  in 
Figs.  48  and  49  show  that  after  the  initial  diffusion  growth  fills 
the  cavity,  the  air  becomes  quite  transparent  and  forces  an  isothermal 
condition  even  through  the  shocks  -  rising  and  falling  uniformly  and 
only  slightly  as  the  shock  strikes  the  wall  and  then  the  origin. 

In  contrast,  the  density  profiles  of  Figs.  50  and  51  show  the 
gradients  in  density  which  travel  with  the  isothermal  shock,  indicating 
the  persistence  of  the  reverberations  even  under  these  radiatively 
transparent  conditions. 

The  final  figure  (Fig.  52)  is  a  composite  of  several  curves  which 
describe  the  transition  from  radiative  growth  to  shock  growth  for  the 
fireball  in  the  one-tenth  normal  atmosphere.  The  shock  density  develop¬ 
ment  is  illustrated  by  the  curve  labeled  (  10tj  -  1),  which  is  the  same 
as  the  compression,  or  the  shock  overdensity  divided  by  the  initial 
ambient  density.  This  curve  suggests  that  a  shock  begins  to  form  by 
12  m  (~  30  p,sec)  but  is  not  fully  formed  until  21  m  (250  y,sec) .  The 
growth  here  should  be  compared  with  that  shown  in  Fig.  12  for  normal 
density  air. 
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The  temperature  at  this  same  front  is  also  shown.  It  drops  most 
rapidly,  of  course,  at  the  time  that  the  compression  rises  the  fastest 
The  inner  temperature  of  the  nearly  isothermal  fireball  is  shown  as  a 
function  of  the  time  only.  The  compression  in  the  bomb  or  case  shock 
is  also  shown,  but  not  at  the  same  time  scale  as  used  for  the  air 
shock  compression  curve,  since  at  corresponding  times,  the  bomb  shock 
is  still  inside  and  behind  the  main  shock  until  the  time  of  catch-up 
(before  200  (j,sec)  . 
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Fig.  1— Shock  position  versus  time, ^  = 
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Fig.  2— Wall  pressure  versus  time.  20  m  cavity, 
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Fig.  3— Wall  temperature  versus  time,  20  m  cavity, 


Fig.  4— Wall  density  ratio  versus  time,  20  m  cavity 
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Fig.  5— Pressure  profiles  at  times  before  first  reflection 
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Fig.  6— Pressure  profiles  at  times  after  first  reflection, 
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Fig.  7  — Pressure  profiles  at  times  before  second  wall  reflection, 
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Fig.  8 — Temperature  profiles  at  times  before  first  reflection, 


-20- 


Fig. 9— Temperature  profiles  at  late  times  (after  many  reflections), 


-22- 


n.  .11  ui  ni  .it  i  u,iuiiuiiiuiiiiiuiuiiiinimtiiiiittuiiuitimiiiii(iiiiiiiiii||iiiimittiuiiii<iiii  n 


uiiiiiiiiitiintiii  niitiitiiiiiiiiiiiiiiiiiiiiiiiltiilfllf  litKttifl 


R  (m) 


Fig.  11  —  Late  density  profiles  (after  many  reflections)/5.^ 
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Fig.  12  —  Radiation  diffusion — shock-wave  transition 
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Fig.  13  — Shock  radius  versus  time  for  40  m  cavity,  / 
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Fig.  16— Wall  density  versus  time  for  40  m  cavity,  p.  = 
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Fig.  17 — Pressure  versus  radius  before  and  after  first  reflection,  P.=P 
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Fig.  18— Pressure  profiles  at  times  between  first  and  second  reflections, 
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Fig.  19  Pressure  profiles  at  late  times, 
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Fig.  20— Temperature  profiles  at  times  before  first  reflection, 


Fig. 21— Temperature  profiles  after  first  reflection, 
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Fig.  22— Density  profiles  before  first  reflection,  p.  =  p q 
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Fig.  23— Late  density  profiles, 


Fig.  24— Wall  pressure  versus  time,  20  m  ca 
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Fig.  25— Wall  temperature  versus  time,  20  m  cavity,  p.  =  p  1 10 


-37- 


Fig.  26— Wall  density  ratio  versus  time,  20  m  cavity,  p.  =  p  1 10 
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Fig.  27— Pressure  profiles,  20  m  cavity,  A  =  P  I 10 
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Fig.  28— Late  pressure  profiles,  P{  =  PQI  10 
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Fig.  29  —  tarly  temperature  profiles,/3—^  ,10 
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Fig.  31  —  Early  density  ratio  profiles  P-^  PQI 10 
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Fig.  32  — Late  density  ratio  profiles,  p.=  p  j  10 


Fig.  33 — Shock  radius 


Fig.  34— Wall  pr assure  versus 
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Fig.  35— Wall  temperature  versus  time,  40  m  cavity,  p.  =  p  1 10 
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Fig.  36— Wall  density  ratio  versus  time,  40  m  cavity,  p.  =  p  / 10 
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Fig. 37 —Pressure  versus  radius  at  times  up  to  first  reflection  P.*Pj\Q 
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Fig.  38— Pressure  versus  radius  at  times 
between  first  and  third  reflections,/3. *P  / 10 
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Fig.  39— Temperature  profiles  at  early  times  (before  first  reflection), 
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Fig.  41  —  Density  versus  radius  at  times  up  to  first  reflection, 
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Fig. 43— Wall  pressure  versus  time,  40  m  cavity,  P-{m  Pj  100 


time,  40  m  cavity,  p.  - p  /100 
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Fig. 45— Wall  density  ratio  versus  time,  40  m  cavity,  P{  =  P01 100 


57 


Fig. 46— Pressure  versus  radius  at  times 
before  first  reflection,  40  m  cavity,  p.^pj  100 
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Fig.  47  — Pressure  profiles,  />.  =  pj  100 
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Fig. 48— Temperature  versus  radius  at  times 
before  first  reflection,  p.  •Pj  100 
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Fig. 49— Temperature  profiles, 
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Fig.  50— Early  density  ratio  profiles,  p.  =  p  i  100 
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Fig.  51  — Late  density  ratio  profiles,  ^.  =  ^Q/100 
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Fig.  52— Radiation  diffusion  transition  to  shock  o.  =  p  / 10 
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